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Abstract 
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1. Introduction 

The fluctuation-dissipation theorem (FDT), vahd for dynamic variables in 
equilibrium, is usually written in a form which involves generalized susceptibilities 
and spectral densities, which are frequency-dependent quantities [l]-[4]. However, 
in order to discuss certain time-dependent properties, it may be more convenient 
to have at hand a formulation of the theorem in the time domain. This need is for 
instance well illustrated in the discussion of aging effects in response and/or cor- 
relation functions of out-of-equilibrium dynamic variables. Then, the equilibrium 
FDT is a priori not applicable and has to be modified by the introduction of a 
violation factor or of an effective temperature, which can conveniently be defined 
in terms of time-dependent quantities [5]-[8]. 

Following arguments outlined in [3], we develop below in a direct and simple 
manner diff'erent ways in which the equilibrium quantum FDT can be formulated 
in the time domain. The corresponding expressions of the theorem, which are 
established in the whole range of temperatures, allow in particular for a discussion 
of both the classical limit [3] and the zero-temperature case. In this latter situ- 
ation, an interesting connection, which does not seem to have been put forward 
previously, is shown to exist with analytic signal theory [9]- [10]. 

2. Time domain formulation of tiie equilibrium FDT 

Let us consider here a system in equilibrium described in the absence of ex- 
ternal perturbations by a time-independent hamiltonian Hq. In the following we 
will be concerned with equilibrium average values which we will denote as (...), 
the symbol (. . .) standing for Trpo ■ ■ with po = e~^^° /Ti e~^^° . 

Since we intend to discuss about linear response functions and symmetrized 
equilibrium correlation functions generically denoted as XBA{t,t') and CBA{t,t'), 
we shall assume that the observables of interest A and B do not commute with 
Hq (were it the case, the response function XBA{t, t') would indeed be zero). This 
hypothesis implies in particular that A and B are centered : {A) = 0, {B) = 0. 

Generally speaking, the time domain formulations of the equilibrium FDT 
establish the link between the linear response function XBA{t,t') and the sym- 
metrized equilibrium correlation function CBA{t, t') = ^[{^{t')B{t)) + {B{t)A{t'))] 
(or the derivative dCBA{t,t')/dt'). 

2.1. The response function in terms of the correlation function 

Consider two quantum-mechanical observables A and B with thermal equi- 
librium correlation functions verifying the property 



{Ait'-inp)B{t)) = {B{t)Ait')), 



P=l/kT, 



(2.1) 



and compute the contour integral 




(2.2) 
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where t is a real time and F the closed contour in the complex r-plane represented 
on Fig. 1. One checks easily that the integrand in / does not present any singular- 
ities inside F. Actually, denoting by |A) and Ex the eigenstates and eigenenergies 
of the system hamiltonian Hq, one has 

{A{T)B{t)) = 1 ^ AAvSA'Ae-^^-+^(^-*)(^--^-'), Z = IVe-^^°, (2.3) 

A, A' 

which displays the fact that the continuation to the region —f3 < $5m r < of the 
complex T-plane of the correlation function {A{T)B{t)) is analytic [3], [11]. 

According to the Cauchy theorem, the integral / is thus equal to zero. In 
the limit R — > oo, one obtains, by gathering the various contributions to /, the 
relation 



/oo „ -1 

dt" [{Ait")Bit)) + {Bit)A{t"))] — , (2.4) 

-oo smh ^ ' 

where the symbol vp denotes the Cauchy principal value. Taking into account the 
Kubo formula for the response function, one gets from Eq. (2.4) the expression 

XBA{t. t') = ^ 0(* - t') VP / dt" CBA{t, t") — , (2.5) 

■^/^ J-oo sinh ^ — - 

where 0(t) denotes the unit step-function. 

Eq. (2.5) allows to compute the response function xba in terms of the sym- 
metrized correlation function Cba- Then, introducing the dissipative part ^ba of 
XbAi as defined by 

XBA{t. t') = 2z Q{t - t') isAit, t'), (2.6) 

one gets from Eq. (2.5): 

iniBA{t,t') = - VP / dt"CBA{t,t") — . (2.7) 

J-oo smh ^ ^ 

Since the two-time equilibrium averages involved in ^ba and Cba only depend on 
the time differences involved, Eq. (2.7) can be rewritten as a convolution product, 
that is 

inisAt)='-CsA(t)*,,^^, (2.8) 
the convolution in the r.h.s. being taken with respect to t. 
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2.2. The correlation function in terms of the response function 

The expression of the correlation function in terms of the response function 
can be derived, either in the same way as above (i.e. by using contour integration), 
or by inverting the convolution product (2.8). 

Compute on the contour F (Fig. 1) the integral 

J=l{Airm))^cot^'-^dr. (2.9) 

Using similar arguments as above (i.e. noticing that J = since there are no 
singularities of the integrand of J inside F), one obtains the relation 

l[{A{t')B{t)) + {B{t)A{t'))] = 



27T 

that is 



7 n 7r(f" - f') 

- VP dt" mt)A{t")) - {A{t")Bm — coth (2.10) 



CbaM = - A VP r dt" [xBAitX)-XABit",t)] ^ COth ^^^^^ (2.11) 

27r ' (in fin 



or 



CBA{t,t') = --VP r dt"iniBA{t,t") coth ^'^ . (2.12) 

TT fin fin 

Before going further, let us add a comment. The arguments which have 
previously been used in the computation of the integral / must be refined in order 
to show, first, that the integrals in the r.h.s. of Eqs. (2.10)-(2.12) are properly 
defined, and, second, that the contributions to the integral J (Eq. (2.9)) of the 
two vertical segments of abscissas R and —R of the contour F (Fig. 1) do vanish in 
the limit i? — > oo. The question stems from the fact that the function coth.{7:t / fih) 
tends towards a finite limit for large values of its argument, contrary to the function 
1/ smh.{TTt / fih) involved in the computation of /, which tends towards zero (this 
latter property insuring that the integrals in the r.h.s. of Eqs. (2.4), (2.5) and 
(2.7) are properly defined and that the contribution to / of the above-mentioned 
segments is actually zero in the limit — > oo). 

A sufficient condition is lim.ji_,oo{-Ai^R — 'i'^y)B{t)) = {A){B) (= since A 
and B arc centered)^. Thus, when treating systems with a finite number of degrees 
of freedom, and oscillating response and correlation functions, in which case this 
limit does not even exist, it will be convenient to introduce a small damping 
which will be eventually let equal to zero, or, which amounts to the same, to treat 
the response and correlation functions as distributions. Further details will be 
provided when treating the harmonic oscillator example (Section 2.4). 



^ Note that for = this condition amounts to limfi^oo{A{i:R)B{t)) — while 
ioT y = P it amounts to lim.ji^oo{B{t)A{±R)) = 0. 
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Eq. (2.12) can be viewed as the reciprocal of Eq. (2.7), since it allows to 
compute the symmetrized correlation function Cba in terms of the dissipative 
part ^BA of Xba- It can be rewritten, using convolution product notations, as 

~ 1 ~ TT iri 

CBA{t) = ih^BA{t) * vp coth— . (2.13) 

TT pn pn 

Note that this latter expression can also be obtained directly by inverting the 
convolution product (2.8), which is easily done by using the relation 

Trrvp— — ^ * — vp coth— = -7r^5(t), (2.14) 
(5% sinh|| (3% ^ (3% ^ ^' ^ ' 

demonstrated in Appendix A. 

Eq. (2.8) together with the inverse relation (2.13) constitute a formulation of 
the equilibrium FDT in the time domain. 

2.3. Relation with the usual frequency domain formulation 

Eqs. (2.8) and (2.13) just correspond by Fourier transformation to the usual 
fluctuation-dissipation relations between the dissipative part ^b^(<^) of the suscep- 
tibility and the Fourier transform Cba(<^) of the symmetrized correlation function 
[l]-[4]: 

iBA{io) = ^ tanh CsAiio), Cba(<^) = U coth C-Ba(<^)- (2-15) 

Indeed, taking the Fourier transform of formulas (2.15), with the follow- 
ing definition of the Fourier transformation F{uo) — J^^dt F{t) e"^^^, F{t) = 

2^ Cba{<^) e~*'^*, and making use of the Fourier formulas (^.3) and (^.6), 

one obtains Eqs. (2.8) and (2.13), rewritten below for clarity: 

~ 1~ TT 1~ 1~ TT Trt 

i^^BA{t) = -CBA{t) * Vp , CBA{t) = iJl^BA{t) * ^ Vp COth— . 

TT pn smh ^ TT pn pn 

(2.16) 

2.4. A basic example : the harmonic oscillator 

Let us consider as a basic example a harmonic oscillator of mass m and angular 
frequency ujq, in thermal equilibrium at temperature T. The oscillator displace- 
ment being denoted by x, one has 

{x{t)x) = [(1 + n) e-^'^o* + n e^'^°*] (2.17) 

and 

{xx{t)) = \n e-''^^* + {l + n) e''^°*] , (2.18) 
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where n = l/{e'^^^° — 1) is the Bose-Einstein function at temperature T. One 
deduces from Eqs. (2.17) and (2.18) the expressions of ^xx{t) and Cxx{t) '■ 



h 



sin (jjQt 



and 



^xxV') = 7^ COth — - — COSCJoI- 



2mujQ 2 

The fluctuation-dissipation relations (2.8) and (2.13) respectively read: 

1 (5?IUJQ TT 1 

smcuot = - COth cos cot * ^ vp 



(2.19) 
(2.20) 



(2.21) 



and 



f3fUjjQ 1 . TT ^ TTt 

coth — - — coso^o^ = sma^o^ * -T^r vp coth — . 

2 TT pn pn 



(2.22) 



Both formulas can be easily checked by a direct calculation, which we report 
in Appendix B, together with supplementary details^ on the contour integration 
needed to compute the integral J (Eq. (2.9)). 

Now, coming back to the general case, let us discuss, first the classical limit, 
then the zero-temperature case. 

3. The classical limit 

Before entering the discussion of this limit, we shall propose another form 
of the first FDT relation (Eqs. (2.7) or (2.8)) in the time domain. It gives the 
expression of the linear response function in terms of the derivative of the equi- 
librium correlation function [3] , and reveals to be especially useful when studying 
the classical limit. 

3.1. Expression of the response function in terms of the derivative 
of the correlation function 

Integrating by parts, one can recast Eq. (2.7) into the equivalent form 



•ft^ f,,'^ 1 r ..n dCBA{t,t") , 
lh^BA{t,t) = - dt — log 



coth 



TTit" - t' 



2(3n 



(3.1) 



or, making use of convolution product notations since ^ba and Cba only depend 
on the time differences involved. 



~ 1 dCBAjt) . 

iTi^BA{t) = * log 

TT at 



coth 



TTt 



2pn 



(3.2) 



^ See the remarks in Section 2.2. 
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an expression equivalent to Eq. (2.8). Thus, at any temperature, the dissipative 
part of the response function (i.e. the function ^BAit)) appears to be proportional 
to the convolution product taken with respect to t of the functions dCsAit) /dt 
and log I coth(7rt/2/3^) |. This latter function is very peaked around t = at high 
temperature while it becomes more and more spread around this value as the 
temperature decreases. 

3.2. The classical limit 

In the classical limit, making use of the property (A. 14), that is 

log I coth ax\ ~|a|^oo 70*^*^^^' ^^-^"^ 
with a = n/2ph, one shows that Eq. (3.2) reduces to 

in is A w = -^ — — ' (3-4) 

which yields for the associated response function XBA{t,t') the expression 

XBA{t,t') = pQ{t-t')^^^g^. (3.5) 

Eq. (3.5) constitutes the expression of the classical equilibrium fluctuation- 
dissipation theorem in the time domain [3] . 

3.3. Illustration : the harmonic oscillator 

One easily checks that Eq. (3.4) is actually verified by the functions ^xx{t) 
as given by Eq. (2.19) and Cxx{t) as given by the classical limit of Eq. (2.20), 
namely 

kT 

Cxx{t) = o cos Wot. (3.6) 

mujQ 

4. The zero temperature case 

4.1. The zero-temperature fluctuation-dissipation theorem 

As the temperature decreases, the function log | coth(7rt/2/3?i)| becomes more 
and more spread around t = 0. At T = 0, coming back to the formulation (2.16) 
of the FDT, one gets 

1 ~ 1 ~ 1 ~ 1 

*^CB^(t) = -CB^(t) * Vp-, CBA{t) = *^^BA(t) * Vp-, (4.1) 

TT t TT t 



that is 



11"^- 1 
-in ^BA {t) = - VP / Cba {t') --- dt\ (4.2) 
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and 

CsAit) = — VP / {-inCBA{t'))—-dt'. (4.3) 

Interestingly enough, these Hilbert transformation relations, which constitute 
the formulation of the fluctuation-dissipation theorem in the time domain at zero 
temperature, are formally similar to the usual Kramers-Kronig relations between 
the real and imaginary parts of the generalized susceptibility, except for the ev- 
ident fact that they hold in the time domain and not in the frequency domain. 
Otherwise stated, at T = 0, the quantities CBA{t) and —ih^BA{t) must constitute 
respectively the real and imaginary parts of an analytic signal ZBA{t) with only 
positive frequency Fourier components [9]- [10]. 

4.2. The zero-temperature analytic signal 

Following these lines, consider the signal ZsAit) = {B{t)A). One has 

ZBAit) = CbaH) + nlBA{t), (4.4) 

which displays the fact that ZsAit) has for real part CBAit) and for imaginary 
part —ifi^BAit)- By Fourier transformation, one gets 

Zba {oj) = Cba {oj) + n ^BA (uj) . (4.5) 
At T = 0, the fluctuation-dissipation relations (2.15) reduce to 

^BAiuj) = ^ sgn(a;) Cba{uj), (4.6) 

with 

sgn(a;) = |_^; ^^0; (4.7) 
so that one gets, as expected, 

Zfl^M = {f (4.8) 

Actually, at T = 0, the function ZBA{t) — {B{t)A) has only positive frequency 
components and thus possesses the characteristics of an analytic signal [9]- [10]. 
This implies that the integral definition Zba{'t) — Zba{'^)&~^^'^ dw 127: can 
then be extended in the whole lower half of the complex r-plane"^ (i.e. '^mr < 0). 
In the same way, consider the function YsAit) = {^B(t)). One has: 

YBAit) = CBAit)-n^BAit). (4.9) 

^ This is in accordance with the above noted fact that at finite temperature 
the prolongation to the region —(3 < Q^mr < of the complex r-plane of the 
correlation function {B{t)A) is analytic. 
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At T = 0, making use of the fluctuation-dissipation relations (2.15), one gets 

Thus, the function YsAit) = {^^{t)) possesses only negative frequency compo- 
nents. The integral definition Yba{t) = Yb a {^) e~ ^'^ du /27r can be extended 
in the whole upper half of the complex r-plane (i.e. Q^mr > 0). 

4.3. Other representations of the analytic signal 

Let us here focus on the analytic signal ZsAit) of positive frequency compo- 
nents (similar considerations can be made for YsAit))- 

For Q^mr < 0, one can write, taking advantage of Eqs. (4.6) and (4.8), 

ZBA{r)^^[ du2hCBA{u)e-''^\ (4.11) 

This yields the following representation of Zba{t) for ^mr < in terms of the 
dissipative part ^BA{t) of the response function: 

Zba{t) = - / dt' i-ih^BAit')) -. (4.12) 

Note that, when Q^mr = 0, the integral in Eq. (4.12) must be understood as a 
principal value. 

Eq. (4.12) can in turn be used as a definition of Zba{t) in the upper half 
of the complex r-plane (i.e. r > 0), where the integral definition Zba{t) = 
/o°° -^sA(<^)e~*'^'^cia;/27r cannot be used. It verifies the property 

Z^Air) = ZBAirn- (4.13) 



4.4. Illustration: the harmonic oscillator 

Let us consider once again the harmonic oscillator of mass m, angular fre- 
quency ujQ and displacement x, in thermal equilibrium at T = 0. 

At T = 0, the equilibrium correlation functions Z^xit) = {x{t)x) and Y^xit) = 
{xx(t)) are given by 

and 

Yxxit) ^ e^'^'^K (4.15) 

Thus Zxx{t) is a monochromatic signal of angular frequency ujq and of Fourier 
spectrum 

Zxx{<^) = 7r5(a;-a;o), (4.16), 

777^0 
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while, similarly, Y^xit) is a monochromatic signal of angular frequency —ujq and 
of Fourier spectrum 

h 

^xxM = 7tS{uj + ujo)- (4.17) 

mujo 

The Hilbert transformation relations (4.1) between Cxx{t) and —ih^xx{t) 
read: 

1 1 1 1 , , 

sin = — cos Wot * vp - , coswo^ = sina;ot*vp-. (4-18) 

TT t TT t 

The representation (4.12) of Z^xi^) in the complex r-plane reads: 

Q^mr > 0, 

(4.19) 

^mr < 0. 

5. Discussion and conclusion 

Clearly, the time domain formulation of the equilibrium fluctuation-dissipation 
theorem is completely equivalent to the widely used frequency form of the theorem. 
In the classical limit, the time domain formulation establishes a proportionality 
relation between the dissipative part of the response function and the derivative 
of the equilibrium correlation function. At zero temperature, it takes the form 
of Hilbert transformation relations between the dissipative part of the response 
function and the symmetrized equilibrium correlation function. 

The time domain formulation is of considerable help in discussing out-of- 
equilibrium phenomena. A good illustration of that can be found in the discussion 
of aging effects. For instance, in a recent paper [8], we have studied these effects 
as displayed by the correlation function of the displacement x{t) — x{to) oi a free 
quantum Brownian particle with respect to its position at a given time to- Indeed, 
since diffusion is going on, the variable x{t) — x(to) never attains equilibrium. For 
any times t and t' such that to < t < t' , the displacement correlation function 
Cxx{t, t';to) depends both on the time difference t — t' and on the waiting time or 
age tw = t' —Iq. Since the particle displacement cannot be viewed as corresponding 
to a stationary stochastic process, Fourier analysis and the Wiener-Khintchine 
theorem cannot be used in order to compute Cxx{t-,t';tQ). This quantity must 
thus be obtained through a double time integration of the velocity correlation 
function Ci;v(ti,t2) (which only depends on the time difference ti — t2 since the 
Brownian particle velocity thermalizes and does not age). 

In this situation, for any temperature of the thermal bath, one can write 
a modified FDT relating the displacement response function Xxx to the partial 
derivative dCxx{t,t']to)/dt' , this latter quantity taking into account even those 
fluctuations of the particle displacement which take place during the waiting time 
(the FDT being valid with no modiflcations only when t„ = 0, i.e. when one wants 



Zxx{t) 



n 

h 



-IU)0T 



11 



to relate Xxx and dCxx{ti't''i'tQ)/dt'\^ _^,)- This is rendered possible through the 
introduction of an associated effective inverse temperature /?eff. (or, equivalently, 
of a violation factor X) depending on both time arguments t — t' and [8]. 

The above considerations, which solely rely on the consideration of time- 
dependent quantities, can be extended to other out-of-equilibrium dynamic vari- 
ables of dissipative systems, classical or quantal [5]- [6]. 
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Appendix A : some useful Fourier transforms and convolution 
relations 

A.l. Fourier transform of tanh fBhuj / 2 

Let us set 

Ji = / auj smut tann . (^-l) 
Jo 2 

Using the expansion 

.00 ^ 



with a; = {tt/ (3?i)x, one gets 



h = ^f2 e-(2fe-i)-t//3ft, t > 0, {A.3) 



k=l 

which yields the final result, valid whatever the sign of t: 

TV 1 



One thus has: 



0, t = 0, 

^ da; e-^'^* tanh ^ = -'i I^. [A.b] 
27r ./_^ 2 TT 



A. 2. Fourier transform of coth/3^a;/2 

Similarly, let us set 



±2= I cujJ svnut coth ■ (^-6) 



Using the expansion 

TTX 2 4a; 



, ..„ _ Ax \ - 1 / , „N 

coth — = \ > — — -, (A.l) 



with Lo defined as above, one gets 

T^ 2 



A;=l 
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which yields the final result, valid whatever the sign of t: 



One thus has: 



2itJ_. 



TV , 

0, t = 0, 



diue-'^' coth^ = -^/2. (A.IO) 

2 TT 



A.3. 

From the relation 

tanh -— coth = 1, (^-H) 

one deduces, by Fourier transformation and use of the above Fourier formulas 
(Eqs. (^.5) and (^.10)), the convolution relation 



At T = 0, Eq. (^.12) reduces to 



vp- * vp- = -7r^5(t). (A13) 



A. 4. Some useful relations 

One has 



log I coth ax\ ~|aHoo Tpi 5(2;) (^-14) 



and 



2 

a TT^ 



sinhax ' ' 2\a\ 



5'{x). (A15) 



Formula (A. 14) can be demonstrated in a standard fashion by considering 
log I coth (XX- 1 BiS Bii distribution. Indeed, for any well-behaved function (/>(x), the in- 
tegral log I cothaxl (l){x) dx tends towards (7r^/4|a|) (/)(0) in the limit \a\ — > 00. 

Formulas (A. 14) and (A. 15) can be used to study the classical limit of the 
convolution relation (A. 12). Indeed, in the classical limit, using Eq. (A. 15) with 
a = 7r/2(3h, the l.h.s. of Eq. {A.12) is seen to reduce to -{n'^/2) 6' (t) * sgn{t) , that 
is, to — 7r^(5(t), as it should. 
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Appendix B : the harmonic oscillator case 

B.l. The first fluctuation-dissipation relation 

The first fluctuation-dissipation relation (2.21) expresses for tlie liarmonic 
oscillator the displacement response function in terms of the displacement corre- 
lation function. The convolution product in the r.h.s. can be easily calculated - 
and Eq. (2.22) easily checked - by making use of the expansion 



1 



sinh|| 



fc=i 



-{2k-l)-n:t/l3h 



t > 0, 



{B.l) 



and computing separately each convolution product of the series. 

In spite of their oscillating character, it is not necessary here to treat the 
oscillator displacement response and correlation functions as distributions. Cor- 
respondingly, due to the fact that the function 1/ sinh(7r/;//3^) involved in the 
computation of the integral / (Eq. (2.2)) tends towards zero for large values of its 
argument, the contour integration on the contour F described in Section 2.1 can 
be carried out without problems, even in the limit i? — > oo. 

B.2. The second fluctuation-dissipation relation 

The second fluctuation-dissipation relation (2.22) expresses for the harmonic 
oscillator the displacement correlation function in terms of the displacement re- 
sponse function. In order to compute the convolution product in the r.h.s., it is 
convenient to write 



coth ^ = sgn(t) + (coth - sgn(t)). 



{B.2) 



with 

sgn(t) 

and to make use of the expansion 



+1, t>0, 
-1, t<0, 



{B.3) 



^ nt 
coth — - sgn(t) 



fc=l 

oo 



fc=l 



(BA) 



The flrst contribution to the r.h.s., namely the convolution product 



1 . TT , . 

--smcvot* — sgn(t). 



(S.5) 
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only makes sense as a relation between distributions^. One has: 

1 2 

sina;ot * sgn(t) = — — cosa^o*- (-^-6) 

TT pncoo 

As for the second contribution to the r.h.s., namely the convolution product 

1 TT 7vt 

- - sin ujot * — (coth — - sgn(t)) , (5.7) 

it can be easily calculated by making use of the expansion (-B.4) and computing 
separately each convolution product of the series. One thus gets: 

-- sma;ot * ^(coth - - sgn(t)) = (coth ^ - ^) ^osu.t. (B.S) 

Gathering together results (-B.6) and (-B.8), one obtains Eq. (2.22). Note that, 
interestingly enough, it is only in the computation of the first contribution to the 
r.h.s. of Eq. (2.22), that is, of the convolution product (-B.5), that the considera- 
tion of the oscillator displacement response and correlation functions as distribu- 
tions is required. 

Coming back to the computation of the contour integral J (Eq. 2.9)), one 
writes, for finite R, 

J = 2mC^a:{t - t') + Ja + Jh. (S.9) 

In Eq. (i?.9), 2inCxx{t~t') is the contribution of the two small semicircles centered 
in t' and t' — ihp, and Ja denotes the contribution to J of the two vertical segments 
of abscissas R and —R (Fig. 1). As for J5, it is defined by 

Jb = VP dt" [{x{t")x{t)) - {x{t)x{n)] ^ coth (B.IO) 

J-R PlT' Pll' 

Interestingly enough, while J a and J^, considered separately, do not tend towards 
a limit when i? — > 00 but display instead an oscillatory behaviour, the sum Ja + Jb 
possesses a well defined limit. Indeed one has, for R finite but such that R » t\ ph, 



2m 



Ja — — _ 2/3 COSCJo-R COSCJot, (-^-H) 



^ This amounts to say that the oscillating functions sin ujot and cos c^o^ have to 
be defined through the usual limiting procedures, namely 

sina;o^= lim sina;o^e~^'*', cosa;ot= lim cosa;ote~^'*'. 
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and 

Jfc ~ TT-T COSUJqR COSUqI COSUjQit — t') 

£sin.„(t-0^(coth=<^^-sgn(t"-t')). (B.12) 
When i? — > oo, the sum + Jb possesses a well-defined limit, namely 



2z7r , ,, 

Ja + Jb = T75 COS a;o(t - t ) 



VP r sina;o(t - t") ^(coth ""^^''^ ^'^ - sgn{t" - t')). (S.13) 

Using then formula (-B.5), and the fact that J = 0, together with the expres- 
sion (2.20) for Cxxii)i one checks again, as expected, the fluctuation-dissipation 
relation (2.22). 



17 



Figure caption 

Fig. 1 

Integration contour for the calculation of the integrals / (Eq. (2.2)) and J 
(Eq. (2.9)). 
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